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Abstract. The neutron charge form factor GEn{(l) is determined from 
^^ \ an analysis of the deuteron quadrupole form factor Fc2{q) data. Re- 

cent calculations, based on a variety of different model interactions and 
<*" • currents, indicate that the contributions associated with the uncertain 

^-^ • two-body operators of shorter range are relatively small for Fc2 (q) , even 

^D . at large momentum transfer q. Hence, GEn{(i) can be extracted from 

L-^ \ Fc2 (q) at large q^ without undue systematic uncertainties from theory. 

o 

^ ■ Introduction. Knowledge of the neutron charge form factor Geh is of great importance 

T^ . for an understanding of its internal structure. It is also crucial for the calculation of 

^ ! nuclear charge form factors, since the latter depend on both Geu. and the proton charge 
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form factor Gep- 

Unfortunately, Geh is still rather poorly known. The difficulties encountered in mea- 
suring Geu are twofold: since there are no free neutron targets, Geu has to be measured 
using composite systems, and this leads to complications due to the presence of other nu- 
cleons. In addition, the electron-neutron cross section is dominated by the contribution 
from the magnetic form factor Gmu, thus making a determination of Geh very difficult. 

The traditional approach to determine Geu uses the deuteron structure function A{q), 
to which the deuteron magnetic form factor, and therefore Gmu, contributes negligibly. 
After removing, via theoretical calculations, the effect of the deuteron structure and the 
contributions to the scattering process from two-body currents, subtraction of the Gep 
contribution then allows one to extract Geu- This procedure is sensitive to systematic 
errors in the theory used, particularly those associated with the modeling of shorter-range 
two-body currents, which are still not very well controlled. 

As a consequence, the resulting values for Geu [Hj HI have fairly large uncertainties, 
and are limited to momentum transfers below q^=16 fm~^. This poses serious problems 
for the calculation of form factors of light nuclei, which one would want to calculate for 
the region covered by data, a region that extends to g^=30-100 fm~^. To the extent 



that current parameterizations of nucleon form factors provide sensible extrapolations for 
Geu at large g^, one must conclude that the contribution of Geu, which seems to fall off 
with increasing q^ much more slowly than Gep, becomes very important at these large 
momentum transfers. 

More recently, the exploitation of a new technique to determine Geu has become 
practical: when performing an (e, e'n) coincidence experiment using polarized electrons 
and when measuring the polarization of the target nucleus or recoil neutron, it becomes 
possible to measure an interference term GeuGmti- This approach removes the difficulty 
associated with GMn-dominance, and is much less dependent on the nuclear structure 
of the target nucleus (deuteron or '^He). Several experiments of this type have been 
performed recently ^^lllll- The resulting values for Geu still have relatively large errors; 
they are, however, mainly statistical and thus can be reduced in the future using better 
technology. The limit in q^ is presently 17 fm~^. Two experiments are under way at 
JLAB to extend the g^-range |T^, |TB[. Within the error bars the available results from 
the double-polarization experiments agree with the values determined from the deuteron 
structure function A{q). 

Exploitation of the quadrupole form factor. In this note, we again use elastic 
electron-deuteron data to determine G_Bn(g) at high momentum transfer. The novel aspect 
of the present approach consists in exploiting the quadrupole form factor Fc2{q) rather 
than the combination of monopole and quadrupole form factors represented by A{q), as 
done in the past. 

When using elastic e-d scattering, two sources of theoretical uncertainty must be 
considered, due to the model for the NN-interaction and the contribution of two-body 
currents. We address the two-body currents first. 

Calculations of Fc2{<i) based on a variety of model interactions and currents indicate 
that contributions from two-body currents are relatively small, even at the high mo- 
mentum transfers of interest here. This is consistent with the naive expectation that, 
since -Fc2(?) involves an integral of the product of deuteron S- and D-wave components 
with the spherical Bessel function J2{qr/2), it is presumably less sensitive to two-body 
currents, at least the short-range ones associated with vector-meson exchanges and/or 
transition mechanisms such as, for example, the pn'f operator, whose contributions are 
quantitatively rather uncertain. This is illustrated in Fig. |l| where we show separately the 
contribution associated with the vr-exchange two-body charge operator, as well as that 
including, in addition, the p-meson and pTTj charge operators. The p-meson and pTTj 
contributions have opposite sign, and tend to cancel each other. As a result, the total 
two-body contribution to Fc2{(l) is dominated, up to g^ ~ 40 fm~^, by the long-range 
vr-exchange operator. 

In this context, it is worth noting that, while modern realistic interactions are essen- 
tially phase-equivalent — they all fit the Nijmegen data-base with a x^ per datum close to 
one — they do differ in the treatment of non-localities. Some of them, like the Argonne fig 



model |T^, are local (in LSJ channels), while some others, like the CD-Bonn model |]15| 
have strong non-localities. In particular, the CD-Bonn interaction has a non-local one 



pion-exchange (OPE) component. However, it has been known for some time |T6|, and 
recently re-emphasized by Forest |jT^, that the local and non-local OPE interactions are 
related to each other via a unitary transformation. Therefore, the differences between 
local and non-local OPE cannot be of any consequence for the prediction of observables, 
such as the deuteron electromagnetic form factors under consideration here, provided, 
of course, that two-body currents generated by the unitary transformation are also in- 
cluded. This fact has been demonstrated JTSf in a calculation of the deuteron structure 
function A{q) and tensor observable T2o{q), based on the local Argonne fig and non-local 
CD-Bonn models and associated (unitarily consistent) electromagnetic currents. The re- 
maining small differences between the calculated A{q) and T2o{q) are due to the additional 
short-range non- localities present in the CD-Bonn. The upshot is that, provided that con- 
sistent calculations — in the sense above — are performed, present "realistic" interactions 
will lead to similar predictions for deuteron electromagnetic observables, at least to the 
extent that these are influenced predominantly by the OPE component. This is especially 
true for the Fc2 form factor for which the vr-exchange contributions dominate. 

Because of these considerations, the theoretical uncertainties for Fc2{<l) are small 
(smaller than for v4(g)), which allows us to to determine Geh with smaller systematic 
errors and extend our knowledge of it to larger q. The use of Fc2{<l) has now become 
possible with the measurements of the polarization observable T2Q{q) in electron-deuteron 
scattering. With T20 known up to g^=40 fm~^, the quadrupole form factor Fc2{<i) can 
experimentally be determined up to that g-value. 
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Figure 1: Effect of the vr-exchange two-body charge operator (dashed) and that obtained 
by including the remaining, shorter range, two-body contributions (solid). 



Experimental Fc2{(l)- 
data on electron-deuteron elastic scattering 



In order to determine Fc2{<i)i we have analyzed the world 

- [^. Some 340 data points on e-d 
scattering are available for momentum transfers below 65 fm~^. The cross sections and 



polarization observables are fitted with flexible paranieterizations for Fcoiq), -Fmi(q') and 
-^02(0') I^H- The statistical errors of the data are calculated using the error matrix. The 



systematic errors, which in general are the largest ones by far, have been evaluated by 
changing each individual data set by the quoted error, and re-fitting the complete data 
set. The changes due to systematic errors of the different, independent, sets of data are 
evaluated separately, and added quadratically. The resulting Fc2{,<i) is used below. 

Determination of Geti- In order to extract GEn we compare to the predictions for 
Fc2{<i) from a number of theoretical calculations. These calculations all use NN potentials 
that provide reasonably good fits to the modern scattering data base, and consistent 
two-body currents. We employ calculations using the Paris and Bonn-B potential by the 



Hannover group [0, ^ , the calculation of Forest and Schiavilla [^ based on the Argonne 
V18 NN potential, and the results obtained recently by the Mainz group |^^ using the 
Bonn OBEPQ-B potential. We also employ the results of the calculation by Van Orden, 
Gross and Devine ES] who use an OBE interaction directly fit to the NN scattering data. 



While the first three calculations are based on an essentially non-relativistic framework 
(with relativistic corrections), the calculation [^ starts from a system of coupled nucleon 
and meson fields and, by means of the Foldy-Wouthuysen transformation, derives the non- 
relativistic limit including all the leading order relativistic contributions. The calculation 
of Van Orden et al. starts from the Bethe-Salpeter equation, which has been reduced to 
a quasi-potential equation by assuming that one of the nucleons is on mass shell. This 
calculation is Lorentz covariant and gauge invariant. All calculations include the relevant 
two-body terms. 

In general, these calculations have used proton form factors as given by the Hoehler 
parameterization [^], which explains well the e-p scattering data up to the q^ of inter- 



est here, including the recent G^jp/GA/p-data [0, |5T|, |5^, Q. The calculations of refs. 
p7| , ^] have been carried out using the dipole form factor for the proton, which only 
roughly reproduces the proton data; here we use the calculation of Arenhoevel et al. p^ 



performed with the Hoehler form factors, while the calculation of van Orden et al. has 
been renormalized to the Hoehler proton form factor. All calculations use the Galster [Q 
neutron charge form factor, or the one by Hoehler, which is very close in the range of q^ 
of interest. 

In fig. 0we show the ratio of these theoretical Fc2{q) form factors to the experimental 
ones. This figure shows that for the C2 form factor the different theoretical predictions 
are quite close. The effect of Geu is appreciable at the higher momentum transfers, large 
enough to be extracted despite the differences between the theoretical predictions. 

In order to determine Geui we use the following approach: As the "theoretical pre- 
diction" we use the average of the five calculations discussed above. For the "theoretical 
error bar" we take the quadratically added deviation of the individual calculations from 
the average. The deviation of this average from experiment we then take as an indication 
that the Galster (or Hoehler) Geti used in the calculation is not quite the correct one, and 
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Figure 2: The ratios of theoretical to experimental C2 form factors as a function of 
momentum transfer. For the Paris potential, we also give the ratio obtained setting 



we determine Geu to get perfect agreement between experiment and the theory average. 
The resulting values of Geu, together with the error bars that include both the spread of 
the theoretical predictions anc? the experimental uncertainty on Fc2, is shown in fig. ^. 

Figure ^ shows that the form factors extracted from the C2 deuteron structure function 
are reasonably accurate in comparison with the results obtained from double-polarization 
measurements, and they agree with them in the g^-region of overlap. In comparison to 
the mean values of Geu determined by Platchkov et al. from the deuteron A{q) structure 
function, the present results are somewhat higher in the region above g^=8 fm~^, but 
compatible with them given the spread of the theoretical predictions available to Platchkov 
et al. at the time. The Geu extracted from the C2-data have larger uncertainties at low 
g^, where the CO multipolarity dominates the cross section and where the available T^q 
data are not very accurate. There, the usage of A{q) leads to superior results. 

The determination of Geti from Fc2 extends to larger momentum transfer than all 
previous determinations, which were limited to q^ ~ 16 fm~^. Somewhat surprisingly, 
the extrapolation of the Galster parameterization beyond its limit of validity (g^=16 
fm~^) does quite well in reproducing the data. As pointed out above, double-polarization 
experiments presently under way at JLAB are expected to provide data in this higher-g^ 
region. 

Conclusions. In this note, we have determined the neutron charge form factor Geu 
starting from the data on electron-deuteron elastic scattering. Contrary to previous anal- 
yses, we use the deuteron quadrupole form factor, which is less sensitive to the short-range 
two-body currents that are not well under control. We employ a representative selection 
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Figure 3: The Geu extracted from the C2 data (o). Also shown are the values obtained 
from double-polarization experiments, and the Galster parameterization with its extrap- 
olation into the region not covered by previous experiments (dotted). 

of both non-relativistic and relativistic theoretical calculations to predict deuteron struc- 
ture functions and contributions of two-body currents, thus allowing use to produce a fair 
estimate of the theoretical uncertainties involved in our procedure. Using this approach, 
we for the first time provide data (other than upper limits [pH]) for Geu at large q^. 
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